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Abstract 

We consider the Casimir energy due to a massless scalar field in a geometry of an infinite wedge 
closed by a Dirichlet circular cylinder, where the wedge is formed by 5-function potentials, so-called 
semitransparent boundaries. A finite expression for the Casimir energy corresponding to the arc and 
the presence of both semitransparent potentials is obtained, from which the torque on the sidewalls 
can be derived. The most interesting part of the calculation is the nontrivial nature of the angular 
mode functions. Numerical results are obtained which are closely analogous to those recently found 
for a magneto dielectric wedge, with the same speed of light on both sides of the wedge boundaries. 
Alternative methods are developed for annular regions with radial semitransparent potentials, 
based on reduced Green's functions for the angular dependence, which allows calculations using 
the multiple-scattering formalism. Numerical results corresponding to the torque on the radial 
plates are likewise computed, which generalize those for the wedge geometry. Generally useful 
formulas for calculating Casimir energies in separable geometries are derived. 
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I. INTRODUCTION 



The Casimir effect which was originally conceived as the attraction between parallel 
perfectly conducting plates, may be regarded as due to the fluctuations of the electromag- 
netic field in the quantum vacuum. In the past six decades, this phenomenon has been 
generalized to many different types of fields and to a variety of geometries and topologies. 
Recent reviews of the Casimir effect include Refs. {2, 3, 4, ^, (|. 

In this paper we will illustrate some new features that arise, for example, in cylindrical 
geometries in which the cylindrical symmetry is broken. In the past three decades there have 
been many works on problems possessing cylindrical symmetry, starting with the calculation 



of the Casimir energy of an infinitely long perfectly conducting cylindrical shell 



7j . The more 



physical but a. 
more recently 



so si gnifi cant 



10 



11 



y more involved case of a dielectric cylinder was considered 



12 



Particularly germane to the present work is the 
calculation of the Casimir effect for a scalar field interior and exterior to a cylindrical 5- 
function potential, a so-called semitransparent cylinder [l^]; in the weak-coupling limit, both 
the semitransparent cylinder and the dielectric cylinder have vanishing Casimir energy. 
The infinite wedge is closely related to the cylindrical geometry. This problem was first 



considered in the late seventies 



16 



171 ] as part of the still ongoing debate about how to 



interpret various divergences in quantum field theory with sharp boundaries and whether 
self-energies of objects have any physical significance. Since then, variations on this idea of 



the wedge have been treated by several authors 18, 



19 



20 



21 



221 ] . and reviewed in Ref. [23 ]. 



A wedge with a coaxial cylindrical shell was considered by Nesterenko et al. 24, 



corresponding local stresses were investigated by Saharian and collatorators 



301 ] . The interaction of an atom with a wedge was studied in Refs. 



31 



32 



26 



33 



25| . and the 



34 



28 



351 ] : this 



geometry is that of the experiment by Sukenik et al. carried out more than 15 years ago 36]. 

n 

Recently, Brevik et al. [37J] calculated the Casimir energy of a magnetodielectric cylinder 
intercut by a perfectly reflecting wedge filled with magnetodielectric material. In all of these 
studies the assumption was made that the wedge be bounded by perfectly conducting walls. 

Although wedges defined by perfect conductors or Dirichlet boundaries break cylindrical 
symmetry, they do so in an easily understood way. When cylindrical symmetry is present, 
the azimuthal quantum number v ranges from — 00 to 00 by integer steps. With a perfect 
conductor, which forces the tangential electric field to vanish on the surface, v takes on values 



2 




FIG. 1: A Dirichlet cylinder intersecting with a coaxial wedge made of semitransparent plates. 



which are related to the opening angle a of the wedge, v = urn/a, where m is an integer. 
But what if the wedge boundaries are not perfect? Recently, Ellingsen et al. j^sl considered 
just such a case, where the wedge was defined by the interface between two magnetodielectric 
media, where the interior sector of the wedge had permittivity E\ and permeability /ii, while 
the exterior sector had permittivity e 2 and permeability [i<i- The geometry was completed 
by inserting a perfectly conducting circular cylinder of radius a centered on the wedge axis. 
To assure a finite result, as well as separability of the problem, the further assumption was 
made there that the speed of light in both media was the same: S\^i = £2/^2- In this case the 
azimuthal quantum number had to be determined by a transcendental equation, which was 
implemented in the calculation through use of the so-called argument principle 39], which 
is just the residue theorem. 

In this paper, we further illustrate this nontrivial azimuthal dependence by considering 
a similar wedge geometry, in which the infinite wedge is formed by two planar ^-function 
potentials, making a dihedral angle a € [0, 7r], closed by a coaxial Dirichlet circular cylinder 
of radius a. See Fig. [TJ We calculate the Casimir energy of a massless scalar field subject 
to these boundary conditions. Because that energy is divergent, we compute the energy 
relative to that when the radius of the cylinder is infinite, and when neither or only one of 
the wedge boundaries is present. Thus, we are computing the energy of mutual interaction 



3 



between the three boundaries. The results, which are rather easily found numerically, are 
very similar to those found for the electromagnetic field in a perfectly conducting cylinder 



with a magnetodielectric wedge, as considered in Ref. [38(] . We describe the geometry in 
terms of cylindrical coordinates, (p, 9, z), with the origin lying along the cylinder axis. 

Because the interest in this problem largely lies in the angular dependence, it is natural 
to approach the problem in an unconventional way, in which the reduced Green's function 
refers to the azimuthal, not the radial coordinate. Technically, that approach requires con- 
sideration of an annular region, which we describe in Sec. IIVI This approach also allows 
use of the multiple-scattering technique, and should have application to more complicated 
geometries, such as the interaction between hyperboloids. We can think of the radial planes 
between the concentric cylinders as forming an annular piston, and we have computed nu- 
merically the Casimir attractive torque between those planes. An alternative approach to 
the determination of the Casimir energy for any such angular potential is described in Sec. PVT . 
The radial functions encountered in these wedge problems are modified Bessel functions of 
imaginary order; since these are rather infrequently described in the literature, we collect 
some relevant properties in Appendix [Bl Required integrals over the squared radial functions 
may be evaluated using identities described in Appendix |A] 

II. SEMITRANSPARENT WEDGE 

In this paper we consider a massless scalar model, in which the wedge is described by a 
5-function potential, V(p, 9) = v(9)/p 2 , 

v {9) = X l S{6 - a/2) + \ 2 5{6 + a/2). (2.1) 

This has the diaphanous property of preserving the speed of light both within and outside 
the wedge. This wedge is superimposed on a coaxial circular cylindrical shell, of radius a, on 
which the scalar field <fi vanishes. To calculate the Casimir energy, we can use the formula 



,2 



2fl 

E = 

where £f is the Green's function for the situation under consideration, satisfying 



E = 1 r ^ / (dr) 2u 2 &(r, r; u), (2.2) 



(-V 2 + V(p, 9) - uj 2 ) &(t, r'; u) = 5(r - r'). (2.3) 
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We can solve this cylindrical problem in terms of the two-dimensional Green's function 



G. 

which satisfies 

1 d d 



/°° rlk 



2 _J__^_ v{9) 

pdp^dp K p 2 d6 2 p 2 



G(p,6;p',9') = -5(p-p')5(6-9'), 
P 



(2.4) 



(2.5) 



where n 2 = k 2 — uj 2 . This separates into two equations, one for the angular eigenfunction 

e„(0), 



d 1 



+ v(9) 



e v (o) = u 2 e u (e), (2.6) 

where we have assumed that the azimuthal eigenfunctions are normalized according to 

d0e„(0)e:,(0) = <w; (2.7) 



orthogonality of the eigenfunctions follows from the Sturm-Liouville nature of the problem. 
Now the two-dimensional Green's function can be constructed as 



G( P , e- P ', e 1 ) = e,(e)e:(e')g v ( P , P ') 



(2.8) 



The radial reduced Green's function satisfies 



1 d d 



p op op pr 



9v(p,p) = -S(p-p'). 
P 



The latter, for a Dirichlet circle at p = a, has the familiar solution, 
9u(p,p) = hinp^K^Kpy) - I u (kp)I u (kp) R " [ '"' ' 



g„(p,p) = Iy^p^K^npy) - K u (kp)K u (kp 



I v {ko) 
„ Una) 



p,p <a, 
p, p' > a. 



(2.9) 

(2.10a) 
(2.10b) 



K v {Ka) 

The azimuthal eigenvalue v is determined by Eq. (12.61) . For the wedge 5-function potential 
(12. ip it is easy to determine v by writing the solutions to Eq. (12. 6p as linear combinations 
of e ±lve , with different coefficients in the sectors \8\ < a/2 and ir > \9\ > a/2. Continuity of 
the function, and discontinuity of its derivative, are imposed at the wedge boundaries. The 
four simultaneous linear homogeneous equations have a solution only if the secular equation 
is satisfied: 

4z/ 2 



= D(u) = sm z v(a - 



7T 



1 — r— r— 1 sin z nv — I y~ + y~ ) sin 2txv. (2.11) 



Because we recognize that the reflection coefficient for a single ^-function interface is = 
(1 + 2iv/\ i )~ 1 , implying that 

Re (rrV) = 1 - — -, Im (r^r, 1 ) = — + — , (2.12) 

we see that this dispersion relation coincides with that found in Ref. 38] when the reflection 
coefficient is purely real. Note that the v — root of Eq. f)2.1ip is spurious and must 
be excluded; unlike for the magnetodielectric wedge, there are no v = modes for the 
semitransparent wedge. 

Now using the general formula (12.21) we compute the Casimir energy per length from 

„ 1 f 00 du 2 f°° dk f 00 , , . . , 

Note that we do not need to know the eigenfunctions 0^, only the eigenvalues v. 

The apparent difficulty, that the eigenvalue condition for v cannot be explicitly solved, 



may be resolved through enforcing the eigenvalue condition by the argument principle 39, 



40 



41 



42 



43 



44j . which gives a sum over non-explicit eigenvalues in terms of a contour 



integral around the real line, 



y^=Ll iv (j^ D{ - v ^) m - (2i4 » 

The contour of integration 7 is illustrated in Fig. [21 Thus, we have the expression after 
making the Euclidean rotation, u i(, and converting to polar coordinates, 

C = kcos<£>, k = Ksimp, (2-15) 



d 



s _ 2/ J dKK 6 J dis\^—\nD(v)j J dppg u (p,p). (2.16) 

This formal expression is rather evidently divergent. We are seeking the mutual inter- 
action energy due to the three boundaries, the two sides of the wedge and the circular arc. 
Therefore, we first must subtract off the free radial Green's function without the circle at 
p = a, which then implies 

/ dppg v {p,p) -»• — ]n[I v (Ka)K v (Ka)]. (2.17) 
J Ik ana 

(The familiar form of this expression is quite general, as illustrated in Appendix |A]) We 

further want to remove the term present without the wedge potential: 

D(u) - DM = ^-^P-. (2.18) 
4z/ z sin ix v 

6 




I 

FIG. 2: Contour of integration 7 for the argument principle (|2. 14[) . Shown also are singularities of 
the integrand (|2.16p along the real and imaginary v axes. 



The resulting energy is still not finite. The reason is that it contains the self-energy of a 
single 5-function potential crossed by the circular cylinder. 

Therefore, we still must remove that part of D due to a single potential, which may be 
obtained by setting A 2 (or Ai) equal to zero: 

Di(i/) = 1 - y cot vir, (2.19) 

so the final form of the dispersion function is obtained by replacing 

DM _ DM bM 



D1HD2W 

A1A2 sin 2 v(a — 7r)/ sin 2 vn + Au 2 — A1A2 — 2z/(Ai + A2) cot vtx 
{2u — Ai cot vk) {2v — A2 cot vit) 

(2.20) 

(Although the spurious v = root is still present, it may be checked that this gives rise to 
an irrelevant divergent constant in the energy.) 

It is now easy to see that the integrand in the expression for the energy falls off exponen- 
tially fast for large v in the right-half complex v plane, except along the real v axis, where 
an exponential convergence factor may be inserted. In particular, for v — irj, rj ^> 1, D(irj) 
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FIG. 3: Contour of integration for the r\ integral in Eq. (|2.24p . 



differs only exponentially from unity: 

D(iv) ~ 1 - 7 ^ rre~ 2m . (2.21) 

Then it is permissible to unfold 7 and convert the contour to one running parallel to the 
imaginary axis as shown in Fig. [31 For imaginary u, v — irj, the dispersion functions become 

D x (it]) = 1 + — coth t]tt, (2.22) 
2rj 



and 



£)/-\_ -A1A2 sinh 2 ?7(a - 7r)/sinh 2 ^7r + 4^ 2 + AiA 2 + 2?7(Ai + A 2 ) cothr]7r (2 23) 
^ ~~ (2r] + Ai coth r]7i) (2t] + A 2 coth r]n) ' ^ ' ' 



Because of Eq. ( 12.211) . the resulting expression for the Casimir energy is manifestly con- 

A 



vergent. This can be further simplified by noting that 4- hi D(irj) is odd, which eliminates 



the K v in Eq. (12.17j) . and then yields the expression 



S 8vr 2 a 2 



1 f°° , 9 f°° , ( d , - , \ d K iri {x) 
— dxx 2 drii — In D(irj) — arctan T v ) ' 2.24 
2 a 2 J J \dr) / dx L ir) {x) 



in terms of 



M*) = 77" [W + , (2.25b) 

Z Sill 7T/i 
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where both Li V (x) and Ki n (x) are real for real r\ and x, and 

I iv (x) = ^^[L iv (x) - iK in {x)\. (2.26) 

We should further note that the arctangent appearing in Eq. (12.241) is not the principal 
value, but rather the smooth function in which the phase is accumulated. (Some properties 
of the modified Bessel functions of imaginary order are collected in Appendix IB1) 
Now we turn to the numerical evaluation of this expression. 



III. NUMERICAL EVALUATION OF CASIMIR ENERGY FOR THE SEMI- 
TRANSPARENT WEDGE 

It is actually quite easy to evaluate Eq. (I2.24p . because the J^logl) function is strongly 
peaked for small 77, except for extremely small values of the dihedral angle a. The diffi- 
culty numerically is that Ki V (x) / Li V (x) is an extremely oscillatory function of x for x < 77, 
becoming infinitely oscillatory as x — > 0. For x > rj, the ratio of modified Bessel func- 
tions of imaginary order monotonically and exponentially approaches zero. (Fo r incomplete 



45|, 



46j; see also 



asymptotic information about Bessel functions of imaginary order see Refs. 
Appendix iBl) The function 

h(v) = [ dx x 2 — arctan l — - , (3.1) 
Jo dx L^(x) 

however, is very smooth. (It vanishes at 77 = 0, so the spurious zero mode should not 

contribute.) To evaluate the double integral, we compute h at a finite number of discrete 

points, form a spline approximation which is indistinguishable from h, and then evaluate 

the function 

f°° d 

e(a)= dr]h{rj)— In D {if}), (3.2) 
Jo V 1 — 1 



numerically. (This strategy is similar to that employed in Ref. 381].) The integrand here 

is quite strongly peaked in a neighborhood of the origin. The Casimir energy, with the 
indicated subtractions, is 

* = ^T-2< a )- ( 3 - 3 ) 
The results found by this strategy are shown in Figs. H] and \5 
These graphs are very reminiscent of those found in Ref. 



for the magnetodielectric 



wedge. In particular, we note that the energies are finite for all a, but as A — ► 00, the 



1 2 3 4 5 6 

a 



FIG. 4: Casimir energy for a semitransparent wedge embedded in a Dirichlet cylinder, as a function 
of the dihedral angle a. Shown in order from highest to lowest are the energies (13.31) for Ai = A2 = 
0.5 to 4.0, by steps of 0.5. 
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FIG. 5: Casimir energy for a semitransparent wedge embedded in a Dirichlet cylinder, as a function 
of the dihedral angle a. Shown in order from highest to lowest are the energies (|3.3p for Ai = 1 
and A2 = 0.1 to 2.1, by steps of 0.5. 




limit corresponding to a Dirichlet boundary, the energy diverges as a 
phenomena was observed in Ref. 



or 2tt: the same 



381 ] for the perfectly conducting wedge limit, treated previ- 



ously in Ref. [37|. This energy should be observable as a torque on the two semitransparent 
plates, r(a) = —^-<£(a), which is, as expect 
the apex of the wedge has been subtracted.} 



plates, r(a) = —^-<£(a), which is, as expected, attractive. (The divergence associated with 
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FIG. 6: The annular geometry considered. 



IV. ALTERNATIVE CALCULATION OF CASIMIR ENERGY FOR SEMI- 
TRANSPARENT WEDGE IN AN ANNULUS 

We start from the formula (12. 2p for the Casimir energy in terms of the Green's function, 

E = - I ^2u 2 Tr(&-^), (4.1) 
2i J 2n 

where the trace denotes the integration over spatial coordinates, and we have again sub- 
tracted out the vacuum contribution. The Green's function 5f(r, r') will satisfy the equation 

[-V 2 - uo 2 + V{r)] &{r, r') = 8{r - r'), (4.2) 

while the free Green's function satisfies the same equation with V(r) = 0. Once again 
we specialize to the cylindrical geometry, but now defined in an annulus. Specifically, we 
require that the boundary conditions on the Green's function are that it vanishes at p = a 
and p = b with b > a, that is, it satisfies Dirichlet boundary conditions on two concentric 
circles. (We will see the necessity for both an inner and an outer boundary in the following.) 
If the potential has the form V(r) = v(8)/p 2 then we can use separation of variables to write 
the Green's function as, in terms of the separation constant 77, 

/°° rfh 
R ri (p; lj, k z )R n (p'; W) k z )g,(6, 9'). (4.3) 
-00 ^ 

The geometry we are considering is illustrated in Fig. EJ Note that instead of expanding in 
eigenfunctions of 9 as in Eq. (12.81) . we have expanded in terms of radial eigenfunctions. (This 
cannot be done without the inner boundary — that is, this alternative separation works for 
an annulus but not for a disk.) The R functions are normalized radial eigenfunctions of the 

11 



eigenvalue problem, 



^ ^ ( 2 i 2\ 2 



R v (p; lu, k z ) = r] 2 R v (p; u, k z ), 



(4.4) 



with boundary values R v (a] cu, k z ) = R v (b; u>, k z ) = 0. The g n is the reduced Green's function 
in the azimuthal coordinates that satisfies the equation 

d 2 



d6 2 



+ r} 2 + v(6) 



g v {6,e') = 6{e-6'), 



(4.5) 



with periodic boundary conditions. Finally inserting Eq. ( 14.31) into Eq. ( 14. II) we get an 
expression for the vacuum energy 

E= ^r t 2uj2 rp f dz ^ f pdp R >^^) I™ m,o)-9?{w)]- 

(4.6) 

We can simplify the result considerably if we now make a Euclidean rotation from u — > 
make the substitution ( 2 + k 2 = k 2 , and integrate out the angle in the (-k z plane to get the 
expression for the energy per unit length 



S = 

4:71 



k dn 



v 



d P R 2 v (p;K) / de[ gri (e,e)-g^(e,9)]. 



(4.7) 



A. The Radial Eigenvalue Problem 



We see that we need an expression for the radial integral 



pdpR 2 (p; k) 



(4.8) 



where the R v s are the normalized eigenfunctions obeying the differential equation (I4.4p . The 
normalization is 

[ b ^R 2 (p-K) = l. (4.9) 

Ja P 

To evaluate this integral we will use the identity (1A9|) . The boundary conditions are 
that -R^(a; k) = R v {b; k) = 0; this is only possible for discrete values of r], namely, this is an 
eigenvalue condition for r\. Let R v (r; k) be a solution to Eq. (14.41) which satisfies R v (a; k) = 
for all f] and k. The normalized solution can then be written as 



(4.10) 



12 



where 



N 2 = I ^R 2 v ( P ,k). (4.11) 



a 



Now writing Eq. (14.41) as 

r^4 +KW G)] w=0 ' (412) 

we can see that ^R^(p; k) is a total derivative given by Eq. (1A9j) . (We replace k — > r\ there.) 
The integral ( 14. lip is now trivial to carry out. We see that the value at the lower limit of 
integration is zero by our boundary condition that R v (a; k) = 0, and the second term on the 
right in Eq. (1A9j) at the upper limit is zero by the eigenvalue condition Rr,(b; n) = 0. This 
gives the normalization constant as 

N2 = ^b^ {b]K) ^ {b]K) - (413) 
Now by considering k rather than 77 as the parameter in Eq. (I4.12p we also have from 
Eq. (1A9j) that the desired integral (14.81) is a total derivative, 

J pdpR 2 v (p, k) = -——R^fo «)^(6; «). (4.14) 
So the desired integral given by Eq. (14.81) can be concisely written as 



b „ d 



pdpR^p-.) = -^^^\. (4.15) 



B. Argument Principle 



Now we again use the argument principle (I2.14p , which we previously used for the angular 
eigenvalues; in this case the sum is over the radial eigenvalues, and the eigenvalue condition 
is given by D(rj) = R v {b; k). So we have occuring in the energy (14.71) the form 

4" J a " 2m J ' R(b-K) I K-i-Rnib-K) 



„ -a - .. Rr,(b;n) \ Kjj-R n (b-K) 

1 f r] d 



f In R v (b; K ). (4.16) 

2717 7 7 KOK 

The expression for the Casimir energy per length (14.71) is then given by 

£ = ^h[ ^ dK J vdv(^nRr,(b-, K )] J d9(g v (0,9)-g<®(9,9)). (4.17) 



7 



13 



C. The Radial Solutions 



The differential equation (I4.12p is the modified Bessel differential equation, of imaginary 
order. We need two independent solutions of this equation, which we could take to be 
K iv (Kp) and L iv (Kp), given by Eq. (I2.25p . Now we want to find the solution R v (p; k) that is 
zero for p = a for all values of r\ and k. An obvious solution is 



where 



Rri(p; K ) = K irj (Ka)I irj {Kp) - I iv (na)K iv (Kp) = R- V {p, k) 



~ 1 , T T . sin vix 

h = -{Iu + I-u) = — L,. 
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(4.18) 



(4.19) 



is the function initially called L v in Ref. 38]; here this is a more convenient choice, in that 
both K v and l v are even in v. 



D. Reduced Green's Function 

We also need the reduced Green's function in the angular coordinates. The free Green's 
function is easily found to be 



g f\e,d') = — (- sinh ti\9 - e'\ + C08i * w cosh n\e - e'\ 

1 2?7 \ ' sinh?77r 



(4.20) 



If we assume a single 5-function potential v (8) = X5(8 — a) then the Green's function is 



2rj 



sinh 7] 1 9 — 9' \ + 



1 



2rj sinh r]7i + A cosh rjTi 



2r] cosh r]7T cosh rj \ 9 — 9' \ 



A 



-{cosh 77(27r + 2a-0-0') - cosh 2r/7r cosh r/|^ - ) J, (4.21) 



2 sinh r/n 

which is defined for 9 and 9' in the interval [a, 2tx + a]. 
The quantity of interest, tx(g — g^), is then 

— A(sinh 777rcoshr]7r + 7]7i) 



d9[ 9ri (9,9)-gM(9,9)] 



2r] 2 sinh 7771" (2?? sinh t/tt + A cosh r)7r) 
and this expression can be seen to be a total derivative 

1 d 



(4.22) 



2n+a 



d9 [g v (9, 9) - gV>(e, 9)] = ±.± In ^1 + A cothr/vr 



2r] drj 



In (l + A^ ) (a, a)), 
(4.23) 

which agrees with the result stated in Eq. (I2.22I) . It is precisely of the expected form (IA11I) . 
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FIG. 7: Large radius limit of annular geometry. For large annular radii a and b, with b — a fixed, 
the annular boundaries become indistinguishable from parallel planes. 



E. Casimir Energy 



The final form for the Casimir energy for a single radial 5-function potential in the annular 
region is 

1 



S 



f ( d 
, K 2 dn / dri ( — — In 
16tt 2 z Jo J 1 \8k 



K iv (Ka)I iv (Kb) - Ii V (m)K iv (Kb) 



I 9 i 
x I — In 

or] 



A 



1 H cothr/7r 

2?7 



(4.24) 



This result may also be obtained by the multiple scattering formalism 47], which says 
that 

E = -^TrlnGG^- 1 = — Tr ln(l + G^V), (4.25) 

the latter form being a useful form for a single potential. We see that Eq. f)4.24p exactly 
corresponds to this if we integrate by parts on 77 and k: 



I f°° f f d r 

dnK J drj ( — In K in (Ka)I iri (Kb) — I ir] (Ka)K iri (Kb) J In [l + A</°'(a, a)] . 



(4.26) 



We can check this result by taking the limit as a and b get very large, but with fixed 
distance between the circles. In this limit this result should reproduce the case of a single 
semitransparent plane between two parallel Dirichlet planes, as illustrated in Fig. [71 The 
energy in that case should be 



in 



K 



lull — 

^rjdfj I 277 



(4.27) 
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where rj 2 = k 2 + (rwr/ (b — a)) 2 . If we use the argument principle we can write it as 



g 



16n 2 i J 



K 3 dK / dfj 



d_ 

dfj 



In 



sin \/fj 2 — K 2 (b — 



f] 2 — K 2 



ld_ 

fj dfj 



In 



1 + 



2fj 



. (4.28) 



The square root divided out in the logarithm is present to remove the spurious square-root 
singularity. It should be noted that both of these expressions (I4.24p and (14.281) are divergent, 
but the divergence is simply the self-energy divergence always present with a single plane. 

It is straightforward to prove that Eq. (14.241) reduces to Eq. (14.281) in the appropriate limit. 
The second logarithm in the former becomes, in the limit 77 — > oo, simply ln(l + A/277), so 
that suggests the correspondence 



n 



A 



A 



(4.29) 



And the leading uniform asymptotic expansion of the modified Bessel functions 45|, |46j gives 



V J 



Kb 



V 



Kb 



V 



Kd 



K irj ( 77— ) Ii J V— ) - K iv ( 77— ) Ii V ( V— ) ~ t{z a ) l l 2 t{z b f' 2 sin [ V {f(z a ) - f(z b ))} 



a/2, 



7] J 2f] 



where z a = na/r], z b = nb/r], t(z) = (1 — z 2 ) 1//2 , and 
f(z) = In 



(4.30) 



(4.31) 



t(z) v ' zt{z) 

(The latter function is a continuation of a function usually called rj, but we have already 
used that symbol repeatedly.) This result holds true for z < 1, but an equivalent form, 
obtained by analytic continuation, holds for z > 1. (See Appendix [B]) Then the derivative 
with respect to k term in Eq. (14.241) becomes 



— \n(KI -IK)~— In 

OK OK 



1 



: sin 



(Ka) 



\/ 1] 2 — {ko) 2 



b — a 



so with the substitutions (I4.29[) . and the observation that 

.2 



9 F( V 2 - K 2 ) = -^^-F( V 2 - K 2 ) 



(4.32) 



(4.33) 



drj w k8k 
we exactly recover the Casimir energy for a semitransparent plate between two Dirichlet 
plates K28\\ . 

It is also easy to check that the energy (14.241) agrees with the expression for the energy 
given by the more conventional approach described in Sec. [Til The latter is 



S 



lQn 2 i 



k 2 cLk I dv 



d_ 

dv 



In 



A 

2z7 



COt UTX 



d 



x — — In [I v {Ka)K u {Kb) — I v (Kb)K v (Ka)\ . 
Ok 



(4.34) 
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\ 

FIG. 8: Transformed contour of integration for the v integral in Eq. (|4.34j) . 



This equivalence may be easily shown by seeing that the integrand is odd in u, because 

I u (na)K u (Kb) — I u (nb)K u (Ka) = I u (na)K u (K,b) — I u (Kb)K u (na), (4.35) 

and then rotating the contour 7 from that shown in Fig. [2] to that in Fig. [3j which may then 
be transformed to that shown in Fig. [SJ by changing v to — v on the negative imaginary axis. 
Thus the contour 7 in the v plane is transformed to 7 in the 7] plane appearing in Eq. (14.241) 
(except traversed in the opposite sense), and the equivalence is established. 



F. Interaction Between Two Semitransparent Planes 

If we want to look at an explicitly finite quantity we will need to look at the interaction 
energy between two semitransparent planes. The geometry is illustrated in Fig. [9j We will 
use a slightly different form of the energy for this, based on the multiple-scattering formalism 
43]: 

E = hJ°° ^M 1 "^^^)- ( 4 - 36 ) 

The subscripts on the Vs represent the potentials Vi(r) = XiS(9)/p 2 , and V^(r) = A2<5(# — 
a)/ p 2 . The Green's functions with superscript (i) represent the interaction with only a single 



17 



FIG. 9: Two semitransparent plates in an annulus. 



potential V^. By using Eq. (14.31) . we can greatly simplify the interaction energy to 

i r°° 

^ = ^1 ^^(l-tr^^V). (4.37) 



We already have an expression for given in Eq. (14.211) . Using the latter we can write 



trg^ ) v 1 g^ ) v 2 



A1A2 cosh 2 r](n — a) 



(4.38) 



(2rj sinh r/n + Ai cosh r/n) (2r] sinh r/n + A 2 cosh 7771") 
This exactly agrees with Eq. (I2.23p . 

Using the argument principle to replace the sum we then get the Casimir energy of two 
semitransparent plates in a Dirichlet annulus, the immediate generalization of the Casimir 
energy (14.261) for a single plate, 



S 



, . k dn I dr] — — In 

i-nh ./ J \orj 

7 



d 



x In 1 - 



K irj (Ka)I irj (Kb) - l iTI (Ka)K iTI (Kb) 
A1A2 cosh 2 rj(n — a) 



(4.39) 



(2rj sinh r/n + Ai cosh r/ir) (2rj sinh r/n + A2 cosh r/n) 
A limiting case when a, b — > 00, b — a fixed, should be two perpendicular semitransparent 
planes, a distance d apart, sandwiched between Dirichlet planes, similar to the single plate 
situation treated in the subsection above. A similar formula should then be 



k dn I dr] ( — In 



sin 



a/?) 2 — K 2 (b — a) 



'fj 2 — K 2 



In 1 



AiA 2 e~ 2 ^ 



(277 + A!)(277 + A 2 ) 

(4.40) 

As in the case of a single plate, this limiting form is immediately obtained from Eq. (14.391) . 

Finally, we verify that we obtain the expression (I2.24p for the wedge geometry. To do 
this, we must include the modes exterior to the outer cylinder (with the wedge extended to 
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infinity as in Fig. [T]) and subtract the energy present if the outer cylinder were not present. 
This means that the radial dispersion function determining the azimuthal eigenvalues rj 
becomes 

R v (b; k) = (K iv (Ka)I iri (Kb) - U^K^b)) (4-41) 

The extended annular energy is then given by Eq. (I4.39j) with R v (b; n) — > R v (b; k). We now 
can distort the contour 7 to one lying along the imaginary axis as shown in Fig. [3l irj — > u, 
(because the second logarithm in Eq. (14.391) falls off exponentially fast for Re?7 > 0), and 
then using the small argument expansion, for real 

nr(! )"'"'■ /»~^^(f)" M . «-a (4.42) 

This means for small a and real v the first logarithm in Eq. (14.391) is \n.I v (Kb)K u (K,b), which 
is just what was encountered in Eq. ( 12. 17ft . We then fold the v integral to encircle the 
positive real axis as in Fig. [8] and integrate by parts in k and v. In this way the form (I2.16P 
is reproduced (with D — > D), which leads to the final expression (12.241) . 

G. Numerical Evaluation of the Casimir Energy for Two Dirichlet Planes in an 
Annulus 

The Casimir energy in Eq. (I4.37P is a quickly converging function so it should be easy to 
evaluate. However, it can be difficult to evaluate the t] eigenvalues, which become functions 
of the wavenumber k and a natural number m. We can get around this problem, again, by 
exploiting the argument principle in order to get a contour integral in the complex plane, 
as in Eq. (14.391) . We cannot integrate along the real line because of the poles introduced 
when we use the argument principle, and unlike with the wedge we cannot open along the 
imaginary axis, because the integral then becomes divergent. So a simple choice is then to 
let the integration run along the angles of 7r/4 and —ir/A in the complex 77 plane. Identifying 
R n (b, k) from Eq. ( 14. 18f) . and writing tr g^vigj? V2 = A(r]) we have 



* = X r r d J Re^^Re^ + Im^^Im^ ^ / W(^) 

I M \l-ReA(Vl, 



4tt 2 



2 \ R V~iv\ ^ 
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Here we have used the property that R v * = R* and A{rf) = A*{rf). The value of R^ip, k) 
is obtained as a numerical solution to the differential equation. Using this technique we can 
obtain a numerical energy in about 1 cpu-second. The results of this calculation are found 
in Fig. [TDJ 
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FIG. 10: This figure shows the energy per length vs the angle between Dirichlet plates. The energy 
is scaled by the inner radius a. The three lines represent three different ratios of inner to outer 
radius £ = 1.1, - = 2, and - = 5. 

Again we would like to compare to known results, so Fig. [11] is a graph of the ratio of the 
energies of an annular piston to a rectangular piston of similar dimension. The rectangular 
piston is constructed so it has the same finite width b — a as the annular piston, and the 
separation distance is the mean distance between the annular plates, 

d= ^2sin(f). (4.44) 

The results make a certain amount of physical sense. The energy of the annular piston 
is greater than that of the rectangular piston for small separation because the inner edge 
of the annular piston is closer, and will contribute more to the energy. However as the 
annular piston gets further away, the other side of the piston will start to contribute and 
lower the overall energy. In addition we see that the small ratio piston is much closer to the 
rectangular piston for small separations than a larger ratio, and in the plateau region for 
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small separations, E ann / E rect « 1.04 for b/a = 1.1 vs. E ann / E iect 1.23 for fo/a = 2. These 
numbers are quite closely reproduced by the ratio of the proximity force approximate value 
of the energy for tilted plates to the energy for parallel plates (ignoring the sidewalls) for 
small tilt angles, 



PFA 



1 a 2 



ay 



16b 2 V + b) ' 



(4.45) 
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FIG. 11: This figure shows the ratio of the energies of an annular Dirichlet piston to a rectangular 
Dirichlet piston of similar dimension vs average separation distance between the plates. The 
variable x is the separation distance scaled by the finite size of the piston b — a, x = d/(b — a). The 
two lines represent two ratios of inner to outer radius | = 1.1, and | = 2. For the latter case, only 
the region a 6 [0, n] is shown. 



V. THETA DEPENDENT POTENTIALS 

Instead of considering, as usually done, spherically or cylindrically symmetric potentials, 
in this paper we have been examining potentials depending on the angles. In particular, in 
two dimensions, in order for separation of variables to work, we have been considering the 
operator 

n2 1 ... d 2 1 d Id 2 1 ... 
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with v{6) given in Eq. (12. ip . The advantage of this potential is that a closed form for the 
secular equation can be given, see Eq. (12.111) . But what can be said for other potentials v (9)7 
The relevant equations for the two-dimensional Green's function are still Eqs. (I2.5l) - (12.10bl) . 
In particular, for the angular eigenfunctions we still have 

= e"(9) + (v 2 -v(9))e(6), (5.2) 

where v is the separation constant. The separation constant is determined from the boundary 
condition in 6. For v (8) a smooth potential one imposes periodic boundary conditions, and 
this is what we concentrate on for concreteness. 

Note, that the only information from Eq. (15.21) that enters Eq. (12. 9p for the radial reduced 
Green's function is the separation constant v. From Eq. (15. 2p its square can be considered 
the eigenvalue of 

Le = ~w +v{e) (5 - 3) 

with periodic boundary conditions. For a nontrivial potential v (8) no explicit form of v 
will be known. But also in general a transcendental equation determining the eigenvalues 



can be obtained; we follow Ref. 48j. In order to formulate this equation let H{8) be the 
fundamental matrix of Eq. (15. 2p . That is, let u^\9) and u^\6) be two linearly independent 
solutions of Eq. (15. 2p . With Wu\6) = duu\0)/d9, the fundamental matrix is 

where we choose the normalizations such that H(0) = 1-2x2 ■ With these definitions and 
normalizations, the equation determining the eigenvalues reads 

= D(u) = (1 - 4 1) (2tt))(1 - w^\2u)) - 4 2 > (270^(270. (5.5) 

The solutions to this equation have to be used in Eq. (12.91) . The Casimir energy expressions 
(I2.16P and (I2.17P then remain valid, once the appropriate radial reduced Green's function 
is used and D(u) given in Eq. (15. 5p is substituted. Once v{6) is specified this allows, in 
principle, for a numerical evaluation of the Casimir energy when suitable subtractions are 
performed. 
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APPENDIX A: AN INTEGRAL THEOREM 



It may be useful to see explicitly how the trace of the subtracted reduced Green's function 
turns into a derivative of a logarithm, as in Eq. (12.171) . Consider a Green's function g K (x, x') 
for a one-dimensional problem described by the differential equation 



~~l~P( x )~j n 2 r{x) + q{x) + p(x)V(x) 

\JjJU \AjiAj 



g K (x,x') = 5(x - x'), 



where V is a ^-function potential, 



V(x) = XS(x-c). 



(Al) 



(A2) 



The problem is defined on the interval a < c < b, where at the boundaries g K satisfies 
Dirichlet boundary conditions, 



g K (a,x') = g K {b,x') = 0. 



(A3) 



If the potential V = 0, let the corresponding Green's function be denoted by g K 



(0) 



Let us solve this problem in terms of two independent solutions of the homogeneous 
equation 

r d d 1 

——p(x)- K 2 r(x) + q(x) u K (x)=0. (A4) 

Let A K be such a solution that vanishes at the left boundary, A K (a) = 0, and B K be an 
independent solution that vanishes at the right boundary, B K (b) = 0, and let them be 
normalized so that the Wronskian is 

1 



W[A K ,B K ](x) = A K (x)B' K (x) - B K (x)A' K (x) 



p(x) 



(A5) 
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Then the "free" Green's function is 

g^(x,x') = A K (x<) J B K (x>), 
and the full Green's function has the form 

aA K (x)A K (x'), a < x, x' < c, 
f3B K (x)B K (x'), c < x, x' < b. 
Now, it is easy to prove that 



g K (x,x') = gf ) {x,x') + 



\Bl(c) 



a 



P 



(A6) 



(A7) 



(A8) 



XA K {c)B K {c) + 1 XA K {c)B K (c) + 1 

It is immediate that any two solutions of the differential equation (1A4p u K and w n satisfy 
d_ 

dx 

and therefore the following indefinite integral follows: 



p{x) ^Tx u ^ x) Tk Wk{x) - u ^ x) Tk1Tx w ^ x] 



2Kr(x)u K (x)w K (x), 



(A9) 



dxr(x)u 2 .(x) 



p(x) 
~2^~ 



-u 



( x ) u > ( x ) — \ n UK ( x ) 
k[X)Uk[X) Ok u'Jx) 



(AlO) 



where u' K {x) = 4-u K {x) 



Now we can evaluate the trace of the interaction part of the Green's function, 
tr (g — g^) = / dxr(x) [g(x,x) — g^°\x,x)~\ 



A 



XA K (c)B K (c) + 1 L 
P(c) XAKc)Bl(c) 
2k XAJc)BJc) + 1 



B 2 K (c) / dxr(x)A 2 K (x)+A 2 K (c) / dxr(x)B 2 K (x) 



A' K (o) d ln A' K (c) B' K {c) d ln B' K (c) 



d 



A K (c)dn A K (c) B K (c)dn B K (c) 
d 



In [1 + XA K (c)B K {c)} = — In [l + Xg^{c, c)] 



(All) 



dn 2 dn 2 
This is the expected expression. As shown in Sec. HVEt this is just the expected multiple 
scattering result. In the Dirichlet limit A — > oo, this agrees with the Bessel function result 
(I2.17p . where p(x) = r(x) = x; although in that case the boundary condition is not Dirichlet 
at the origin, p(0) = 0. 



APPENDIX B: MODIFIED BESSEL FUNCTIONS OF PURE IMAGINARY OR- 
DER 



In this work we encountered the following differential equation 

d d 



x—x— x + ri ib(x 

ox ox ' 



(Bl) 



24 



which is the modified Bessel equation, with the wrong sign for the order parameter rf . The 
solutions are then obviously modified Bessel functions of imaginary order. So we might 
choose as the independent pair of solutions the modified Bessel function of the first kind, 
of positive and negative pure imaginary order Ii V (x) and However the J^'s are not 

numerically satisfactory functions; their values for real x are complex, and the phase is x 
dependent. A standard pair of functions can be defined as 

K„(x) = - I v {x)) , (B2a) 

2 sin V7i 

Lu(x) = — (J„(x) + I_ v {x)) . (B2b) 

Z sin UTT 

The K v (x) is the standard modified Bessel function of the second kind, also called the 
Macdonald function. Both K ir] (x) and L iri (x) are real for real values of r\ and x. For a fixed 
77, both K and L oscillate with relatively constant amplitude for x < rj and they die or grow 
exponent ially for x > rj respectively. The limiting behaviors are given in this appendix; see 
Refs. 45l. 1461] . Although this definition of L u is convenient in Sec. HI1 for the considerations of 
Sec. IIVt the sin utt in Eq. flB2bj) introduces spurious singularities, and it is more convenient 
there to simply use 

I u (x) = \ (I u (x) + I_ v (x)) = S -^L„(x), (B3) 
2 in 



also called L v in Ref. 38| . In the following we will give the behaviors of Ki V and J^. 



1. Small Argument 



For fixed 77 > 0, in the limit as x — ► 



K iv (x) 



I if) 



7T 



77 sinh ?77r 
/ sinh?77r^ ' 



1/2 



x 



sin ( rj In — 



+ (x z 



X 



cos ( rj In — 



where <f>„ is given by 



b v = arg[T(H-M7)]. 



(B4a) 
(B4b) 

(B5) 
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2. Large Argument 



For fixed 77 > and large argument \x\ — > 00 we have 

K iv {x) = (JLy /2 e -*[l + 0( x -i)], \ aT gx\ < (B6a) 
/ 1 \ 1/2 

hr,(x) = I — ) e x [l + ^(x- 1 )] , I argx| < - - 5, (B6b) 



for arbitrary 5 > 0. 



3. Uniform Asymptotic Expansion 

The uniform asymptotic expansions are for both large order and argument. For fixed 
z > we have for the leading behavior 

-»pr/2 / AC \ 1/4 

K iv ( V z) ~ [^ 2 j Ai(-^C), (B7a) 

J ^)-2^(l^J Bi(-,n), (B7b) 

where the Ai(x) and Bi(x) are the Airy functions of the first and second kinds respectively, 
and ( is given by the relation 



/Z7T e ' ' / 7T\ 

K iv {r)z) ~ \j — { i_ z 2y/A cos " 4 J > ( B9a ) 



^C 3/2 = /W , /(*) = In ( 1± ^^) - V / T^. (B8) 

For 2 < 1 we can use the behavior of the Airy functions for large negative argument to 
simplify the expressions, 

~ ~V 2^ (l-^i/4 si11 " 4 J • ( B9b ) 

If we choose the branch cut for equation (1B8|) such that ( is a continuous real function of z, 
then for z > 1 we can simplify the expressions to read 

Ki„(r,z) ~ v /^ (z2 _ 1)1/4 e-»W, (BlOa) 

^ ~ ^(i^ 6 ^* (B10b) 
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where g(z) is the natural extension of f(z) 

g(z) = — arcsec z + V z 2 — 1 . 
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